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We derive an effective coupling parameter for two-dimensional Yukawa systems based on the 
height of the first maximum of the pair distribution function. Two variants - one valid in the high- 
coupling range, the other for arbitrary couplings of the liquid - are derived. Comparison to previous 
approaches to Yukawa coupling parameters shows that the present concept is more general and 
more accurate. Using, in addition, dynamical information contained in the velocity autocorrelation 
function, we outline a reference data method that can be employed as a non-invasive measurement 
scheme of the plasma parameters - the coupling strength and the screening length. This approach 
requires only input from a time-series of configuration snapshots and particle velocities with no 
recourse to additional information about the system. Our results should be directly applicable as a 
simple, yet reliable diagnostic method for a variety of experiments, including dusty plasmas, colloidal 
suspensions and ions in traps, and can be employed to facilitate comparisons between experiments, 
theory and simulations. 
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I. INTRODUCTION 

The Yukawa model p] of interaction between charged 
particles is successfully used to describe a large variety 
of experimental systems. Among these are charged dust 
grains in complex plasma liquids [2}^ and colloidal sus- 
pensions of particles 5 ]. Due to its simplicity and uni- 
versality as well as its success in describing experimen- 
tal phenomena, the Yukawa model has been the subject 
of extensive theoretical investigations. Static properties 
such as the triple point and phase transitions have been 
investigated in three dimensions 6 8] and are of contin- 
uing interest in two dimensions where the question of the 
melting scenario is still a topic of some controversy [9j- 
I12j . Furthermore, dynamic processes of Yukawa systems 
under investigation include transport properties such as 
diffusion [T3J, thermal conductivity and viscosity P^HTS] 
in three dimensions. In two dimensions (2D), anomalous 
behaviour of the transport properties has been a subject 
of detailed investigations [PTH2lTj . The collective excita- 
tion spectra and waves are also a subject of ongoing in- 
terest, in three [21] as well as in two [22] dimensions and 
in quasi-2D systems [33J and systems under the influence 
of a magnetic field [2~Hf2l)] . 

It is desirable to test these and other theoretical results 
with data obtained for example from dusty plasma and 
colloidal suspensions experiments or ions in traps [2"Tll2"8] . 
These systems can be realized in a highly correlated 
yet liquid-like state, where particles are able to over- 
come the energy barriers of their local potential wells 
("cages") with a sufficiently high probability to destroy 
long-range order indicative of solid systems. Using ad- 
vanced particle-recognition schemes, these experimental 
setups deliver reliable phase-space trajectories (position 



and velocities) of the individual particles in the correlated 
system, see, e.g., Ref. [3] and references therein. 

In order to compare theoretical results for Yukawa sys- 
tems to experiments, one first has to make sure that 
the experiments are performed in conditions where the 
Yukawa model is applicable. This requires, in particu- 
lar, that ion streaming and wake effects are negligible 
and the pair interaction potential is close to an isotropic 
one. These conditions are rather well known. For exam- 
ple, the particles should be located sufficiently far from 
the electrodes or from void regions. Once this has been 
established [and throughout this paper this will be as- 
sumed] , it is crucial to infer from the experiment the two 
parameters governing the Yukawa model: the Coulomb 
coupling parameter, T, and the inverse screening length 
k (see Section |n] for definitions). For k = 0, the Yukawa- 
model reduces to the well-known one-component-plasma 
(OCP) model. 

In the case of 2D dusty plasma setups, T and k are 
often obtained by comparison of experimental wave dis- 
persions with theoretical and simulation data [29] 130] . 
Another approach is to excite sound waves by laser ma- 
nipulation and measure the ratio of the sound speeds 
of compressional and shear waves which is a function of 
k |31j or the transverse sound speed Ct alone, which, how- 
ever, is not very sensitive to variation in k (AC t /C t < 0.2 
for < k < 2) [321 [33J . The particle charge can alterna- 
tively be measured by observing the oscillation frequency 
of individual particles in the sheath [34j [35] or by theoret- 
ical models of dust charging [SB] . A parameter combining 
r and k non-uniquely into a single value can be obtained 
by observing the particle diffusion as a function of time 
alone [37] . All these methods have a limited range of 
applicability and they are, in part, not easy to use. 
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It is, therefore, desirable, to investigate alternative or 
complementary methods which rely on easily accessible 
experimental data, which is the goal of this paper. Wc 
first develop an effective coupling parameter T* for 2D 
Yukawa systems which can be used to compare struc- 
turally similar systems, i.e. systems with (almost) iden- 
tical pair distributions, but different values of T and k. 
This enables the comparison of Yukawa systems across 
different values of n, and makes it possible to investigate 
the influence of the interaction range alone on the be- 
haviour of a Yukawa system without simultaneously and 
inadvertently disrupting structural features. 

We then propose a reference data method (RDM) for 
obtaining (T, k) simultaneously from equilibrium particle 
trajectories alone, requiring no further information about 
the system as input. Our approach towards the RDM 
consists of first obtaining reference data by molecular dy- 
namics (MD) simulation and then condensing structural 
and dynamical data into classification numbers. Fitting 
the relations between these classification numbers and 
(r, k), we derive symbolic formulas to obtain the system 
parameters. These formulas for F and k are fairly easy 
to use and provide an accuracy of the order of 10 percent 
which is sufficient for most experimental purposes. The 
RDM requires no external perturbation of the system. 

The remainder of this paper is structured as follows: 
In Secti on |Tl} wc introduce our simulation technique. In 
Section |III[ we derive the effective coupling parameter 
and compare it to other approaches. This coupling pa- 
rameter is used in Section ITVl to derive the RDM scheme. 
We conclude in Section IVl 

II. YUKAWA-MODEL AND SIMULATION 
TECHNIQUE 

The Yukawa model is described by a collection of par- 
ticles whose pairwise interaction is given by the potential 

$(r) = - e~ r/XD . (1) 
r 

Here r is the distance between two particles, Ad is the De- 
bye screening length and q is the charge of the particles. 
In thermal equilibrium, the Yukawa system is character- 
ized by two dimensionless numbers: 

i) the Coulomb coupling parameter, T = q 2 /(aksT), and 

ii) the dimensionless inverse screening length n = a/ Ad, 

where (in 2D) a = [nn] , n is the areal number den- 
sity of the particles, and T is the temperature. Be- 
low we will use the (inverse of the) plasma frequency 

1/2 

uip = [nq 2 1 (sqw)] as the time scale and a as the length 
scale. 

We obtain reference data for the Yukawa model by 
solving the resultant equations of motion for N — 
4080 . . . 16320 particles interacting via the pair poten- 
tial by molecular dynamics simulations. The par- 
ticles are uniform in mass m and charge q and are placed 



inside a rectangular simulation box subject to periodic 
boundary conditions. Forces are calculated by impos- 
ing a K-dependent cut-off radius, when k > 0. For the 
Coulomb case of K = 0, we employ the Ewald summation 
technique |38| . The system is put into equilibrium prior 
to data collection by continuous rescaling of the parti- 
cles velocities to the required temperature. During the 
data collection, the system evolves microcanonically. For 
details of the simulation, see Ref. [33] ■ 

III. EFFECTIVE COUPLING AND PAIR 
DISTRIBUTION 

In this section, we study the pair distribution function 
of a 2D Yukawa system for different values of the cou- 
pling parameter and of the screening length. We discuss 
some previous work which introduced effective coupling 
parameters and propose an improved approach which is 
based on the shape of the pair distribution, more pre- 
cisely, on the height of its first peak. In conclusion, we 
demonstrate that the proposed concept can be extended 
from the strongly coupled fluid regime also to the freezing 
point. 

A. Pair distribution function 

The pair distribution function (PDF) is defined as 




where r t j — |rj — rj\, (.) is a temporal average and the 
primed sum indicates the omission of the term i = j. 
It reflects the probability of finding two particles sep- 
arated by r in the system [39]. The PDF of an ideal 
(non-interacting) system is unity, while correlations due 
to interactions give rise to regions of enhanced [g(r) > 1] 
and reduced [g(r) < 1] values of the PDF. 

Figure [ija) shows the PDF in a Yukawa system (n = 2) 
for different values of T. All PDFs show a correlation hole 
at small r due to particle repulsion and different degrees 
of modulations deriving from particle ordering, which be- 
come more pronounced as the coupling T is increased. In 
the liquid phase, the order in the system is only short- 
ranged and the peaks of the PDF decay quickly with 
increasing r. 

The effect of the screening strength, k, is analyzed 
in Fig. [ijb). With increased screening, the correlation 
hole narrows and the modulation of the PDF becomes 
less pronounced. For reasonably high coupling, Yukawa 
systems with different values of n possess a very similar 
PDF, if T is chosen appropriately for each k, cf. Fig.[ljc). 
This behavior was observed by various authors, see, e.g., 
Ref. [22] and references therein, and suggests that the 
structural properties of a Yukawa liquid do not depend 
separately on T and k but rather on a single parameter 
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FIG. 1: Pair distribution function of a strongly coupled Yukawa liquid, a) at constant screening, k = 2, and different F, b) at 
constant Coulomb coupling, T = 120, and different k and c) for four different combinations V, k which yield almost identical 
PDF, in particular a constant height of the first peak. 



effective coupling parameter which is a suitable combina- 
tion of the two. 



B. Effective coupling parameters 

For the OCP, k = 0, the coupling parameter T is 
equivalent to the ratio of the nominal nearest-neighbour 
interaction energy q 2 /a and the thermal energy ksT. 
Clearly, for k > 0, such an interpretation is not possi- 
ble, since the factor exp(— k) is missing from the former. 
Thus, the Coulomb coupling T has no immediate phys- 
ical meaning for systems with n > 0. The observation 
in Fig. [ljc) that Yukawa systems with different combina- 
tions of (r, k) show similar (static) properties has led to 
efforts to define a universal effective coupling parameter 
T* (r, k) analogous to the OCP coupling parameter. Usu- 
ally, T* is tailored to converge to T as k goes to zero, so 
that for example the freezing of a Yukawa liquid occurs 
at the known freezing point of the OCP, T* = 137. 

Various authors have proposed different functional 
forms of r*(r, k). The intuitive definition 

rj(i» = r e - K (3) 

has been proposed by Ikezi and others, e.g. P0rl4"2"] . How- 
ever, one problem with this approach is that systems with 
the same do not exhibit many similarities and can be 
either in the liquid or in the solid phase, depending on 

k [IU1II3]. 

A different definition given by Vaulina and Khra- 
pak [44], 

T* 2 (r, K ) = cre- dK (i + dK+ {d ^j , (4) 

has been frequently used for three dimensional systems, 
where c = 1 [JSJ |551 |3S] ■ In two dimensions, the factor 
c = 1.5 is sometimes used [47] . The factor d in the above 
equations appears due to differences in the respective sys- 
tem of units and is d = n 1 ^ 2 in 2D and d = (47T/3) 1 / 3 
in 3D. Equation Q can be obtained by assuming that 
the coupling parameter is proportional to l 2 <fr"(a)/2 at a 



characteristic distance I instead of being proportional to 
$(a). This expression is sensitive to the distance fluctu- 
ations of particles around a stationary state (character- 
ized by the interparticle distance a) and should thus be 
more suitable to characterize the system behavior near 
the crystallization point than the absolute value of the 
potential itself. A similar expression was used in Ref. [15] 
for finite three dimensional systems. 

For two dimensional systems, Hartman et al. [10l |22"1 
I49j have defined an effective coupling parameter, 

T* 3 (T, K )=f( K )T ) (5) 

based on a constant amplitude of the first peak of the pair 
distribution function g(r). The scaling function /(k) has 
been found by fitting numerical data to a polynomial in 
k up to fourth order which has the form 

f( K ) = 1 - 0.388k 2 + 0.138k 3 - 0.0138k 4 . 

The validity of the definition of Tg is limited i) to the 
maximum value of K max = 3.0 used to derive f(k) [50j . 
and ii) to a minimum value of Tg 40. Below this latter 
value, the scaling is not exclusively dependent on k and 
thus no universal scaling function /(k) can be given, as 
will be shown below. 

The ratio of the coupling parameter and the coupling 
parameter F c (k) at the crystallization transition temper- 
ature has also been employed to quantify the actual phys- 
ical coupling of the system [5Ur[52] : 

ri(r,K) = ^r. (6) 

This definition correctly captures the main structural 
properties of a Yukawa system but its validity range re- 
mains unclear. 

Our goal in this paper is to find an effective coupling 
parameter T* which is valid (at least) in the entire liquid 
phase. It should correctly reflect the structural proper- 
ties of the system which are fully contained in the pair 
distribution function. For strongly coupled systems, the 
first peak amplitude of the PDF already characterizes 



4 



3.0 
2.8 
2.6 
2.4 
2.2 
2.0 
1.8 
1.6 




10 



100 



1000 



10000 



FIG. 2: The height g m ax of the first peak of the PDF as a 
function of T. 



the complete PDF, except for very high values of r (see 
discussion in Ref. [TO])- For less strongly coupled liquids, 
the size of the correlation hole and the position of the 
subsequent peaks of the PDF differ for systems with dif- 
ferent values of (r, k), even if the first peak height of the 
PDF is identical. Using the first peak height to charac- 
terize the structural properties of the system is, of course, 
only one of several possible choices but offers numerous 
advantages; it is, for example, clearly defined, dimension- 
less and easily accessible in simulations and experiments. 

We, therefore, employ a similar approach as Hart- 
mann et al. by finding a scaling function w which results 
in a constancy of the first peak amplitude of the PDF. 
However, to remove the limitations of their formula we 
allow for a T-dependence of the scaling function, 



r*(r,«) = w(t,k)t, 

and extend our definition to 0.0 < k < 5.0. 
we are interested in the inverse relation 



r(r* 



w*(T* 



(7) 

In addition, 
(8) 



As a first step, we obtain - by first-principle MD sim- 
ulations - the reference data for n = ... 5 in steps of 
Ak = 0.5 and 60-90 values of T per n. From these data, 
we read off the PDF peak height <? ma x as a function of T 
and k. Some of the resulting data are shown in Fig. [2] 
We chose an interval g max G (1.5,3.0) to carry out the 
subsequent analysis. 

The ratio g m ax(r, k) normalized by g m ax(r, 0) is equiv- 
alent to w*(T*,k). Employing symbolic regression [53] . 
we find that the following functional form provides excel- 
lent agreement with the simulation data: 



exp 



A + B(T*)k 
B(T*) = Bx+Bs/r*, 



(9) 
(10) 



where A = 2.37221, B x = 0.55515, and B 2 = 5.56585. 
For reasonably high coupling, T* > 40, the function 
B(T*) is only weakly dependent on T* (cf. Fig. [3]) and 




FIG. 3: B(T*) according to Eq. (10 1. The horizontal line is 



the constant approximation used in Eq. (Ill (see text 



is thus well approximated by a constant. This leads to 
a simpler approximate effective coupling which depends 
on k only, in line with the earlier definitions ([3|-([6]), 



F 



w (k)i , 



r* > 40 



(11) 



where w*(n) = w*(T, k), as defined above but with A = 
2.33425, Bi = 0.63529, and B 2 = 0, resulting from the 
described averaging procedure. 

This relation is readily inverted, yielding 



w(k)T . 



(12) 



where w(n) = 1/w*(k). While the general solution of 
Eqs. ([8]), ([9]) for V* is quite challenging and can only be 
achieved numerically, we have found that the following 
symbolic expression can be used as an approximation, in 
lieu of the numerical inversion: 



w(T,k) = C(k) + 



D(k) 



D{K) = (dJ + d 2 ) 



(13) 
(14) 
(15) 



Here d = 0.0063503, C 2 = 0.00209988,C* 3 = 0.0668851 
and D x = 0.309364, D 2 = 1.749. 



C. Analysis of formula (111 



Let us first compare the r*-independent definition ( fTTj ) 
with the alternative scaling functions which were pro- 
posed earlier, cf. Eqs. ([3])-([5]). They are shown in Fig. [4] 
together with the actual ratio of the peak heights ob- 
tained from the MD simulations which are indicated by 
symbols. Obviously, the simulation data for different val- 
ues of r* do not fall on a universal curve; therefore, the 
scaling cannot be a universal function of k. How well a 
particular scaling function fits the data is dependent on 
the interval of T* which it aims to describe. We also note 
that the exponential and Vaulina's definition are not, un- 
like Hartmann's and ours, specifically aimed at describing 
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FIG. 4: The ratio F/F* as a function of k. Symbols denote 
data taken from Fig. [2] and the full line corr esp onds to the 
scaling function according to the fit formula of this pa- 
per. The expression |5| of Hartmann et al., the formula Q 
of Vaulina et al. (c = 1), and the expression by Ikezi are 
shown for comparison. The thin dotted line is the continua- 
tion of |5| beyond its scope of validity. 



FIG. 6: The ratio F/F* as a function of ft. Symbols denote 
data taken from Fig. [2] and the full lines correspond to the 
scaling function according to our full model Q. 




FIG. 5: K-dependence of the ratios of pairs of different scaling 
functions, (a) Our formula (111 to formula (JsJ [Hartmann et 
al.], (b) Formula |TT| to formula Q [Vaulina et al.], (c) For- 
mula ([5| to formula Q. The ratio (c) continues to decay 
monotonically beyond k = 3; the horizontal lines mark a de- 
viation of ±5%. 



the structural order and are thus not expected to coincide 
with the data points of Fig. |4j 

Apart from the exponential scaling, the curves 
nonetheless show a nearly universal behaviour. Below 
K s» 2.5, there are essentially only small differences. 
Vaulina's scaling is slightly higher than the other two 
and fits the structural data somewhat worse. Differences 
between Hartmann's et al. and our scaling are small and 
likely due to the T* interval of the respective validity, as 
discussed above. Note however, that our scaling oj*(k) 
contains only two numerical fit parameters as opposed to 
Hartmann's three. 

Above k = 3, Hartmann's scaling is not valid any more 
and its numerical extrapolation, shown by the thin dotted 
line in Fig. [4] clearly gives unphysical results. Vaulina's 
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FIG. 7: The ratio of the F*-independent scaling function 
w*(k) to the full expression w*(T*,k) for three values of k. 
The horizontal lines mark a deviation of ±5%. 



and our scaling differ appreciably in this range of k. 
To quantify the respective deviations, various ratios of 
pairs of scaling functions are shown as a function of k 
in Fig. [5] It can be seen that the difference between 
Hartmann's scaling and ours only slightly exceed 5% for 
k < 3 while differences between Vaulina's and our scal- 
ing quickly grow with increase of k and exceed 15% for 
k > 3. 



D. Analysis of formula Q 

We derived the r*-independent scaling by averaging 
the full model expression over its r*-dependence for 
r* > 40. As can be seen in Fig. |4j this procedure gives 
best agreement for T* « 80 but it fails for smaller val- 
ues of r*. Thus, for systems with lower coupling, the full 
scaling function ^ should be used instead. In Fig. [6j this 
scaling function is shown for four values of T* covering 
all but the most extreme states of liquid systems. The 
excellent agreement is evident, even though only three 
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free parameters enter into the scaling function. 

Finally, we directly compare our two models: the ap- 
proximate r* -independent scaling function, Eq. ([Ill, 
and the full T*-dependent counterpart, Eq. ([8]). Figure 7] 
shows, for three values of k, their ratio as a function of T* . 
As expected, deviations are smallest for T* s» 80 . . . 100 
and increase with increasing distance from this range. 
The relative deviations are smaller for small k and are 
less than 3% for k = 1 and T* e [20, 120]. For higher k, 
the deviations increase substantially as the spread of the 
ratio r/r* becomes larger (cf. symbols in Fig. E|. 

The effective coupling parameter introduced in this 
section can serve a number of purposes. It can be used 
to compare systems with different interactions but iden- 
tical structural properties by varying the screening from 
k = (long-ranged interaction) to n — 5 (short-ranged 
interaction). Such a line of research promises to give in- 
sight into the dependence of physical properties on the 
range of interparticle interaction. Structural properties 
such as the crystallization temperature can be obtained 
by taking the well-known results for the OCP and trans- 
ferring them, via T* , to Yukawa systems. Since T(r*,K) 
grows monotonically with k, it is also possible to study 
dynamical properties of systems with identical structure 
but different thermal excitation (recall that T oc T -1 ). 
This allows one to separate the influence of thermal ex- 
citation and structural changes on dynamical properties, 
which arc otherwise changed simultaneously when T is 
varied 126 . 



E. Pair distribution and crystallization 

While our present analysis focuses on strongly cou- 
pled Yukawa fluids it is, nevertheless, tempting to take 
a brief look at larger couplings where the system freezes. 
The change of the pair distribution function (and of re- 
lated functions, such as the static structure factor) at the 
liquid-solid transition has been studied in great detail by 
many authors [54 -57]. Yet we are not aware of a system- 
atic analysis of the height of the first peak of g(r) at the 
phase transition. We have, therefore, performed exten- 
sive additional MD simulations for T and k values in the 
vicinity of the freezing transition. 

Some representative results for the first peak height 
<?max of the PDF are shown in Fig. [8] and extend the 
data of Fig. [2] to larger couplings. The first observation 
is that <7 max changes discontinuously. Comparing to the 
known freezing points of Yukawa and Coulomb systems 
we observe that the position of this jumps exactly agrees 
with these values. We thus conclude that the height of 
the first peak of g{r) may serve as valuable indicator of 
the liquid-solid transition in one-component Yukawa sys- 
tems. This value is fairly easily detected in an experiment 
and does not require long measurements to achieve con- 
vergence. Preliminary results indicate that this behavior 
is observed also for other types of interactions and for 
phase transitions in quantum systems [55H60J as well. 
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FIG. 8: The height <? ma x of the first peak of the PDF as a 
function of T in the vicinity of the liquid-solid transition for 
a Coulomb system and a Yukawa system with n — 1. The 
vertical lines mark the crystallization points, T c — 137 and 
r c = 187, respectively. 



IV. THE REFERENCE DATA METHOD 

In the previous section, we have introduced an effective 
coupling parameter T* which is based on the constancy 
of the first peak of the PDF, i.e. on the universality 
of short-range structural features. However, even if two 
systems are practically indistinguishable in their struc- 
tural properties, and thus have identical T*, in general, 
these systems will exhibit differences in their dynami- 
cal behaviour. We can reformulate this statement in the 
following way: If we could measure structural and dy- 
namical properties, this information should be sufficient 
to distinguish Yukawa systems with respect to both their 
coupling and screening parameter. 

We will now exploit this fact to derive a practical 
method which allows to calculate both the values of T 
and k simultaneously from a time-ordered series of par- 
ticle snapshots alone. No assumptions about the time 
scale of the snapshots are being made, except that they 
are evenly spaced in time (a restraint easily accomplished 
in experiments, but easily dropped if needed) and of suf- 
ficiently rapid succession to resolve single-particle oscilla- 
tions. This allows for computing the instantaneous par- 
ticle velocities. 

As a first step, we obtain a direct mapping between 
the peak height g max of the PDF and the effective cou- 
pling parameter introduced earlier. We find the following 
simple relation: 



r*( ffmax ) - 19.2581 g 2 v 



16.6829 g n 



(16) 



which allows one to calculate T* without prior knowl- 
edge of k and r [and so far without knowledge of the 
dynamical properties] of the system. Simply by calcu- 
lating g(r) in a very limited range of r, as can be done 
from a couple of system snapshots, one can obtain the 
effective coupling which governs the system's static and 
structural behaviour. 

What remains to be achieved is to extract T and k from 
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FIG. 9: Time dependence of the velocity autocorrelation 
function Z{t) for systems with identical short-range order, 
i.e. identical F*, cf. Fig. [T^c). (a): Time argument normal- 
ized to the nominal plasma frequency, (b) : Time normalized 
to the Einstein frequency. Curves in (a) are shifted for clar- 
ity [Z(0) = 1 in all cases]. The arrow in (a) indicates the 
definition of Az (see text). 



the result ( 16 ) which does require additional information. 
A suitable supplementary quantity is the velocity auto- 
correlation function which we discuss in the following. 



A. Velocity autocorrelation function 

The (normalized) velocity autocorrelation function 
(VACF) Z(t) is defined as 



Z(t) = 



(v(t) ■ v(0)) 
(\v(0)\ 2 ) 



(17) 



where (.) is an ensemble average which is performed by 
averaging over all particles and v denotes the particle 
velocity. 

The VACF is a fundamental measure of the dynamics 
of an ensemble of particles and is intimately connected, 
e.g., to the diffusion process. In reasonably highly cou- 
pled liquids, at any time instant, a fraction of particles 
oscillate in their local potential cages before diffusing far- 
ther. These periodic trajectories manifest themselves in 
an oscillatory VACF, which decays with time because of 
anisotropies of the cages and the increasing number of 
particles which leave their potential cage. 

Figure [9]ja) shows the VACF for systems with identi- 
cal effective coupling T* but different re. The oscillation 
period of the VACF curves differs appreciably across dif- 
ferent values of re. However, when replotted as a func- 









■p* on 






X\120 


40 

GO 

sn 








ou 

100 — 




r* = 2iN\ 




120 











0.1 0.2 0.3 0.4 0.5 0.6 
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(Eq. ( 18 1) for different values of T* . 



tion of time t multiplied by the Einstein frequency uje, 
cf. Fig. [9^b), the VACF oscillation frequency becomes 
very similar |10j . This is not surprising, since loe is de- 
fined as the (average) oscillation frequency of a single 
particle in the frozen environment of all other particles 
and thus takes into account the range of the pair inter- 
action. 

Our goal is to use the dynamical information contained 
in the VACF to distinguish between systems with the 
same short-range structural order (i.e., systems with the 
same V*) but different re. If one would have knowledge 
of the nominal plasma frequency lo p of the system, this 
would be easily accomplished by using the re-dependence 
of loe or, equivalently, of the VACF oscillations. How- 
ever, uj p is usually not known and cannot be calculated 
from the trajectory snapshots alone. Therefore, we must 
choose a criterion which is independent of the time scale. 
While the periodicity of the VACF can be mapped on a 
single curve for different K, Fig. [9jb) shows that, never- 
theless, there are clear re-dependent differences in the am- 
plitude of the oscillations. This suggests to use, e.g., the 
amplitude difference Az between the first minimum and 
the second maximum of the VACF as a parameter suit- 
able to distinguish between different screening parame- 
ters, see Fig. [9]ja). 

Using our reference data as before, it is possible to 
derive an explicit symbolic approximation for the screen- 
ing parameter re as a function of the two "experimental" 
parameters Az and r*(g max ): 



A A (n = 
f A (n = 

B A (T*) = 
C A = 



A; 



A z A z - C A 

/A(r') 



B A (T*) 

1.28486 - T 
0.00171645 - 0.011499 F 
0.166741 - 0.0013331 T* 
0.856117. 



(18) 



This result for n(Az,T*) is shown in Fig. 10 as a func- 
tion of Az, for various values of T* . The somewhat con- 



voluted form of Eq. ( 18 1 results from a strong sensitivity 



of the functional form of the curves on their parameters. 





r sim 


K 


r 





60 


-0.09 


65.93 


1 


100 


1.12 


109.0 


2 


300 


2.07 


323.8 


3 


300 


3.01 


297.8 


4 


1000 


4.07 


1062 


5 


3000 


5.06 


3246 



TABLE I: Examples for the application of the RDM. The 
actual parameters of the simulation carry the subscript "sim" 
(first and second column), the results from the RDM are in 
the third and fourth columns 



^sim 


(|Ak|> 


(|Ar/r sim |) 





0.10 


0.12 


1 


0.15 


0.10 


2 


0.09 


0.03 


3 


0.05 


0.12 


0...5 


0.12 


0.12 



TABLE II: Accuracy of the RDM for different values of K s im, 
averaged over V* . The second (third) column contains the 
absolute (relative) deviations of the RDM results from the 
original simulation parameters. 



Figure [TU] also demonstrates that Az is a viable choice to 
characterize Z(t) as a function of T* and k: The curves 
do not cross and proceed monotonically both for fixed T* 
and fixed k. In addition, the values of Az span a suf- 
ficiently broad interval to make an accurate separation 
possible. 



B. Applying the RDM 

The proposed reference data method is now a straight- 
forward application of the results derived before. It can 
be written down as follows: 

1. From the particle snapshots, calculate g(r) and 
Z(t) and read off g max and Az- 



2. Calculate T* using Eq. ( 16 ) and g 



3. Use Eq. (18) to calculate K from T* and Az- 



4. Calculate T using Eqs. ^ and ^ or Eq. ((TTj 



To assess the accuracy of the RDM, we use it to calculate 
r and k from the g max and Az values of our MD simu- 
lations. Some typical examples are shown in Tab. [TJ The 
full results (Tab.|llj) indicate an average absolute error of 
An < 0.15 (recall that k is given in dimensionless units) 
and a relative error Ar/r < 0.15 for the experimentally 
most relevant values of k < 3. 



V. CONCLUSION 

In this work, we have derived an effective coupling pa- 
rameter r* for strongly correlated 2D Yukawa systems 
which is based on the constancy of the first peak of the 
PDF. For high coupling, T* > 40, a simple one-parameter 
relation holds, T*/r = w(k). This expression was tested 
against previous definitions of an effective coupling pa- 
rameter and shown to give similar (k < 3) or even sub- 
stantially improved (k > 3) results in terms of short- 
range structure. 

In the low coupling regime, however, no such simple 
form exists, and we showed how to define an appropriate 



coupling parameter which is of the form r*/T = w(k, T). 
It again allows for an adequate representation of the 
structural features of Yukawa systems at these couplings. 
Only two (three) free parameters enter the proposed 
functional form for high (low) coupling. 

In the second part, we have outlined a reference data 
method to calculate the values of T and n of a Yukawa 
system from a time series of configuration snapshots. Us- 
ing the dynamical information contained in the first os- 
cillation cycle of the VACF, we were able to devise a 
method which does not depend on additional physical in- 
formation such as temperature, mass, charge or density, 
to achieve this goal. We tested this method and found 
that it gives reliable results with an error not exceeding 
10 . . . 15% for majority of relevant situations. It should 
be directly applicable to a number of systems of charged 
particles, in particular dusty plasmas. 

We mention two caveats one should be aware of when 
applying the results of this work: 1) the effective cou- 
pling parameter and the RDM were derived for PDF 
peak heights in the range of <7 max = 1.5 . . . 3.0 which cor- 
responds to r* 18 . . . 124), and are thus restricted to 
this range; 2) the effect of friction on the dynamics of 
the system has not been considered. This is not of im- 
portance for r* but more so for the RDM. Future work 
will aim to relax these restraints. 

Furthermore, we have shown that the peak height of 
g(r) is also sensitive to the liquid-solid transition and 
may serve as a simply accessible first estimate. Finally, 
while our results were derived for spatially homogeneous 
macroscopic 2D Yukawa systems we expect that they can 
be straightforwardly extended to systems in an external 
trap such as atoms, molecules, ions or dusty plasmas. 
This is based on the observation [61] [62] that the shell 
structure of these systems is accurately described within 
a theory which uses the pair correlations of the associ- 
ated homogeneous system. This question is the subject 
of ongoing work. 
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